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Abstract 

The accurate knowledge of luminosity at e^e~ flavour factories re- 
quires the precision calculation of the Bhabha cross section at large scat- 
tering angles. In order to achieve a theoretical accuracy at the 0.1% 
level, the relevant effect of QED radiative corrections is taken into ac- 
count in the framework of the Parton Shower method, which allows ex- 
clusive event generation. On this scheme, a Monte Carlo event generator 
(BABAYAGA) is developed for data analysis. To test the reliability of 
the approach, a benchmark calculation, including exact 0{a) corrections 
and higher-order leading logarithmic contributions, is developed as well, 
implemented in a Monte Carlo integrator (LABSPV) and compared in 
detail with the BABAYAGA predictions. The effect of initial-state and 
final-state radiation, 0{a) next-to-leading and higher-order leading cor- 
rections is investigated and discussed in the presence of realistic event 
selections. The theoretical precision of BABAYAGA is estimated to be at 
the 0.5% level. 



1 Introduction 

The precise determination of the machine luminosity is necessary for the suc- 
cessful accomplishment of the physics program of e"'"e~ colliders operating in 
the region of the low- lying hadronic resonances, such as DA$NE (Frascati) ||l| , 
VEPP-2M (Novosibirsk) g, as weU as for the BABAR [|, | and BELLE |, | 
experiments at PEP-II and KEKB. In particular, the precise measurement of 

iTalk given by CM. Carloni Calamc at the workshop PHYSICS AND DETECTORS FOR 
DA<I>NE (Frascati, 16-19 November 1999), to be published in the workshop proceedings. 
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the hadronic cross section requires a luminosity determination with a total rel- 
ative error better than 1% It is well known that the luminosity of e'^e~ 
colliders can be precisely derived by the relation L = N/ath, where N and ath 
are the number of events and the theoretical cross section of a given reference 
reaction. In order to make the total Imninosity error as small as possible, the 
reference process should be characterized by a large cross section and calculable 
with high theoretical accuracy. At low-energy e'^e~ machines, the best candi- 
date fulfiUing the above criteria is the Bhabha process (e+e^ e+e~) detected 
at large scattering angles. 

On the theoretical side, precision calculations of the large-angle Bhabha 
(LABH) cross-section are therefore demanded, with a theoretical accuracy at the 
O{10~^) level. This requires the inclusion in the calculation of all the relevant 
radiative corrections, in particular the large effects due to photonic radiation. 
The complete and exclusive simulation of events in generators is also strongly 
required by the experimental analysis. 



2 Theoretical approach 

The calculation of the Bhabha scattering cross-section, corrected by the effects 
of photon radiation, and the corresponding event generation is performed ac- 
cording to the master formula 

cr corrected = J dx^dx+dy^dy+ J dQiabD{x^,Q'^)D{x+,Q'^) X 

)e{cuts). (1) 

In the previous equation, the electron Structure Function (SF) D{x,Q^) is the 
solution of DGLAP equation in QED. It takes into account soft-photon expo- 
nentiation and multiple hard bremsstrahlung emission in the leading log (LL) 
approximation both for the QED initial-state (ISR) and final-state radi- 
ation (FSR). The QED-DGLAP equation can be exactly solved by means of 
the QED Parton Shower (PS) algorithm which allows also exclusive pho- 
ton generation in the LL approximation. A more detailed discussion about 
the implementation of the PS algorithm as adopted in the present analysis will 
be given elsewhere |10 . In eq. |l|, dao/dQ is the "hard-scattering" differential 



cross section relevant for centre of mass (cm.) energy around 1 GeV, includ- 
ing the photonic s- and i-channel diagrams, their interference and the (small) 
contribution due to $ exchange. In the hard-scattering cross section, the cor- 
rection due to vacuum polarization is taken into account as well, according to 
the parameterization and the recipe given in ref. [0| . The effect of the running 



2 



coupling constant at y/s ~ is to enhance the cross section by « 2(2.5)% 
for 20°(50°) < i?± < 160°(130°). The jacobian factor j{x-,x+,9iab) in eq. g 
accounts for the boost from the cm. to the laboratory frame due to emission 
by initial state e"*" and of unbalanced radiation, while Q{cuts) stands for 
(arbitrary) experimental cuts implementation. 

Upon the above-sketched theoretical background, a new Monte Carlo (MC) 
generator (BABAYAGA) for simulation of the LABH process at ^-factories 
has been developed. In the program both ISR and FSR are simulated and 
the complete kinematics of the generated events is reconstructed in the LL 
approximation. The possibility of performing an up to 0{a) calculation of eq. |^ 
is included as well, in view of a comparison with the exact 0{a) perturbative 
results. 

In order to test the precision and the reliability of the Bhabha generator, 
an exact 0{a) calculation has also been addressed, by computing the up to 
0(a) corrected cross-section, consisting of soft -I- virtual jl^ and hard photon 
corrections [|3j. Moreover, higher-order LL terms can be summed on top of the 
exact 0{a) cross section whitin the collinear SF approach, following the algo- 
rithm of ref. [Q. This formulation is available in the form of a MC integrator 
(LABSPV), which is a suitable modification of the SABSPV code described in 
ref. 



3 Numerical results 

In the following, the selection criteria adopted for the analysis and the simu- 
lations correspond to realistic data taking at DA<I>NE and VEPP-2M, at cm. 
energy y/s = 1.019 GeV. The energy cut imposed on the final-state electron and 
positron is Ei^^^ = 0.4 GeV, with the angular acceptance of 20° < ■d± < 160° 
or 50° < i9± < 130° and the (maximum) acoUinearity cut allowed to vary in the 
range £,max = 5°-25°. 

A sample of simulations obtained by means of BABAYAGA is shown in fig. 0, 
where the energy, the cosine of the angle, the p± of the most energetic photon 
of each event and the missing mass of the event are plotted. As expected, the 
behaviour of photonic radiation (soft and collinear to charged particles) is well 
reproduced by the PS. The effect of FSR has also been investigated. The cor- 
rected cross section including only ISR has been compared with the corrected 
cross section including both ISR and FSR. We noticed that, as a consequence of 
the rather severe cuts, the total effects of photon radiation is to reduce the inte- 
grated cross section by an 0(10%) amount. As expected, half of the whole effect 
must be ascribed to FSR when non-calorimetric ( "bare" ) event selection fll, f6i 
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Figure 1: Differential distributions obtained by means of BABAYAGA. Clock- 
wise: energy, cosine of the angle, p± of the most energetic photon and missing 
mass. 

is adopted. 

The comparison between the exact 0(a) calculation and the 0(a) predic- 
tions of the PS generator allows to evaluate the size of the 0(a) next-to-leading- 
order (NLO) corrections missing in the LL approximation PS predictions. More- 
over, this comparison can be a useful guideline to improve the agreement be- 
tween perturbative and PS results, for example, by properly choosing the vir- 
tuality in the electron SF in such a way that the bulk of 0(a) NLO terms is 
effectively reabsorbed into the LL contributions. The scale choice = st/u (s, 
t and u are the usual Mandelstam variables) allows to keep under control the 
dominant structure due to initial-, final- and initial-final-state interference radi- 
ation [0. As a function of the acoUinearity cut, the relative difference between 
the exact 0(a) cross section and the corresponding PS one is shown in fig. ||, 
for the angular acceptances 20° < 'd± < 160° and 50° < 'd± < 130° and for two 
different choices of the scale in the PS, i.e. = st/u and = 0.75 • st/u. 
It can be seen that, with the scale = 0.75 • st/u, the difference between the 
exact 0(a) calculation and the PS predictions is within 0.5%. This naive ex- 
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0(a) test - BABAYAGA vs LABSPV 
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Figure 2: Relative differences between the exact 0{a) and the up to 0{a) PS 
cross sections, as functions of the acoUinearity cut and for two choices of the 
scale in the electron SF's, in two different angular acceptance regions. 



ample illustrates how, for a given selection criterion, the level of agreement can 
be substantially improved by a simple redefinition of the maximum virtuality 
of the electromagnetic shower. Going beyond this simple recipe would require a 
true merging between perturbative calculation and PS scheme, which is beyond 
the scope of the present analysis. 

In addition to the evaluation of the 0{a) NLO corrections, it is important, 
for an assessment of the theoretical precision, to quantify the amount of the 
higher-order LL contributions with typical experimental cuts. The size of LL 
0(a"L") {n > 2) corrections can be derived in the PS scheme by comparing the 
full all-order predictions with the corresponding up to 0{a) truncation, as shown 
in fig. ^. The comparison shows that the 0(a"L") corrections are unavoidable 
for a theoretical precision better than 1%, being their contribution 0.7% at 
^max = 5° and 0.3-0.4% for larger acoUinearity cuts in the angular acceptance 
20° < i9± < 160°. The 0{a"L") corrections are even more important, of the 
order of 1.5%, in the narrower angular range 50° < i9± < 130°. It is worth 
noticing that the generators so far used by the experimental groups at Frascati 
and Novosibirsk include only jl^ 0{a) corrections, missing the important effect 
of higher-order contributions. 
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Figure 3: Relative deviation between 0{a) and all- orders corrected cross section, 
in the PS scheme, as a function of acollinearity cut. Two different angular 
acceptances for e+ and are considered. 



4 Conclusions and perspectives 

In order to provide predictions of interest for the luminosity determination at 
e~^e~ flavour factories, a precision calculation of the LABH process has been 
addressed. It is based on a QED PS algorithm (the details of the formulation 
as adopted in the present paper will be given elsewhere |lO[), which accounts 
for corrections due to ISR and FSR (and interference) in the LL approxima- 
tion and allows the complete event generation. A new MC event generator 
(BABAYAGA) has been developed and is available for a full experimental sim- 
ulation; actually, it is under test at Frascati and Novosibirsk. The overall preci- 
sion of the PS approach has been checked by means of a benchmark calculation, 
which includes exact 0{a) and higher-order LL corrections and is available as a 
MC integrator (LABSPV), allowing for precise cross section calculations. Crit- 
ical comparisons between the exact 0{a) and the 0{a) PS calculations pointed 
out that the contribution of the 0{a) NLO corrections is important for the re- 
quired theoretical precision. Moreover, the effect of higher-order 0{a"'L") LL 
corrections has been evaluated to be at the 1-2% level. 
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By virtue of its generality, the PS approach could be employed to simulate 
and to evaluate radiative corrections to other large-angle QED processes, as for 
example e^e~ — > 77 or e^e~ — > fi'^fi^. An interesting application of PS would 
be the simulation of processes with tagged photons, e.g. e+e^ hadrons + 7. 

In conclusion, our analysis points out that theoretical predictions aiming 
at a O(10~^) precision must include the effects of both 0{a) NLO terms and 
0{a"L'^) LL contributions. As a consequence of that, we can estimate the 
present accuracy of our generator BABAYAGA to be at 0.5% level and the 
accuracy of the integrator LABSPV at 0.1% level. In the future, an improvement 
of the presented approach is needed by means of an appropriate merging of the 
exact 0{a) matrix element with the exclusive photon exponentiation realized 
by the PS algorithm. 
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